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ABSTRACT 

We explore correlations between visual extinction and polarization along the 
western side of the IRAS Vela Shell using a published polarimetric catalog of sev- 
eral hundreds of objects. Our extinction maps along this ionization front (I-front) 
find evidence of clumpy structure with typical masses between 1.5 and 6 M Q and 
a mean length scale L ~0.47 pc. The polarimetric data allowed us to investi- 
gate the distribution of the local magnetic field in small (~pc) scales across the 
I-front. Using the dispersion of polarization position angles, we find variations 
in the kinetic-to- magnetic energy density ratio of, at least, one order of mag- 
nitude along the I-front, with the magnetic pressure generally dominating over 
the turbulent motions. These findings suggest that the magnetic component has 
a significant contribution to the dynamical balance of this region. Along the 
I-front, the mean magnetic field projected on the sky is [0.018 ± 0.013] mG. The 
polarization efficiency seems to change along the I-front. We attribute high po- 
larization efficiencies in regions of relatively low extinction to an optimum degree 
of grain alignment. Analysis of the mass-to-magnetic flux ratio shows that this 
quantity is consistent with the subcritical regime (A < 1), showing that mag- 
netic support is indeed important in the region. Our data extend the overall 
A— A r (H 2 ) relation toward lower density values and show that such trend contin- 
ues smoothly toward low iV(H 2 ) values. This provides general support for the 
evolution of initially subcritical clouds to an eventual supercritical stage. 



Subject headings: polarization — ISM: clouds — ISM: individual (Gum Nebula, 
Vela Shell) — ISM: extinction — ISM: magnetic fields 
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Introduction 



Magnetic fields are impor tant for the dynamics of expanding shells (jSpitzerlll978uTroland fc Heiles 



shells with magnetic fields exist in severa 
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(Hartauist. PiliDD. & Havnes 
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L994; 


Redman et al. 1998; 


Wardle 




1998; 


Williams, Dvson & Hartquist 


1 2000l: 


Williams & Dvson 


2001 


) . A proper determination 





of the magnetic field in shells is necessary to constraint these models. 



Following IChandrasekhar fc Fermil (119531 . hereinafter CF), the dispersion of polarization 
position angles, observed through a polarizing medium, permits to estimate the strength of 
the magnetic field component in the plane of the sky (B). In addition, this disper sion also 
yield s information about the kinetic-to-magnetic energy density ratio (pkin/Pmag, 
1990h . 



Zweibel 



Previous works used the CF procedure to calculate the magnetic field projected on the 



sky in several regions of ISM (IGonatas et ah] Il990l ; iMorris et al.l Il992l ; IChrysostomou et al. 



1994 ; lltoh et al.lll999l ; iHenning et al.ll200ll ). Nevertheless, a relatively large number of back 
ground stars are necessary in order to obtain a reasonable estimate for the strength of the 
magnetic field in a small area. Good statistics requires a sampling d ensity higher enough 
to ex ceed that of the large-scale structures of the local magnetic field (lAndersson fc Potter 
2005h . 



Polarimetric catalogues with several hundreds of objects i n selected ISM regions and 
higher sa mpling density are now available for this type of study (jPereyra fc Magalhaesll2002 
- Paper 1. 12004 ) with typical value of ~1.2 objects/arcminute 2 that is, at least, three orders 
of magnitude higher than previous works. 

The IRAS Vela Shell (IVS) encloses a cavity that appears to have been formed by st ars of 
the V ela OB2 association through the effects of stellar wind s and supernova explosions (ISahul 
1992h . However, the existence of the IVS is controversial and Woermann. Gaylard. fc Otrucek 
r l2ooih suggested that it is instead a density enhancement in the Gum nebula. 



A section of the western side of IVS (see Figures 2 and 3 in Paper I) shows a well- 
defined ionization front (I- front) seen almost edge-on with several magnetic field patterns. 
Churchwell et al. Jl996h found that the kinetic energy is one order of magnitude higher than 
the gravitational potential energy in this region. They concluded that the structure is not 
gravitationally bound and would disperse on very short timescales were it not for the ram 
pressure of the expanding IVS, which continually sweeps up new interstellar matter into the 
cloud. 
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Fig. 1. — Extinction map toward the western side of IRAS Vela Shell. The gray contours 
go from 0.3 to 1.0 mag in steps of 0.1 mag. In black contours (from 1.0 to 1.2 mag) are the 
extinction structures (in letters) with Ay > 1. 
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In this work we explore the role of the magnetic field on the dynamics of the IVS region 
using a published polarimetric catalogue (Paper I), with 856 objects. In §2], we describe the 
technique used to construct the extinction maps. We also detect extinction structures and 
estimate their masses. In §31 we present a statistical analysis of the polarization catalogue. 
In §H we compare the polarization along the I-front with the detected extinction structures. 
In §5j we apply the dispersion of polarization angle to infer the distribution of pun/ Pmag 
and the strength of B along the I-front. In §6j we explore the polarizing efficiency and 
its implications about the amplification of the local magnetic field. In §7J we analyze the 
magnetic support in terms of the mass-to-magnetic flux ratio. Conclusions are drawn in §5J 
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Table 1. Extinction structures. 



ID RA 20 oo DEC2000 Ay -/V(H 2 ) area L mass obs. 

( h : m : s ) ( : ' : " ) (mag) (10 21 cm- 2 ) ( /2 ) (pc) (M©) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



A 


07:43:26.90 


-42:22:22.55 


1.42 


1.33 


8 


0.52 


3.73 V12 a , 


B 


07:42:52.65 


-42:24:43.90 


1.16 


1.09 


4 


0.37 


1.68 B b 


C 


07:42:43.71 


-42:16:16.90 


1.05 


0.99 


4 


0.37 


1.52 


D 


07:40:05.85 


-42:58:11.84 


1.50 


1.41 


12 


0.64 


5.95 



A 1 



Vilas-Boas. Mvers. & Fuller! (11994! ) 



IChurchwell et all jl99fih 
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2. Visual Extinction by automatic star counts 



2.1. The extinction map 

We used our own tasks to obtain maps of visual extinction of the region using 

an automatically star count technique. Firstly, we obtained positions of objects with stellar 
profiles from a Digitized Sky Survey image (see Figure 3 in Paper I) using the IRAF task 
daofind. This routine automatically detects objects over certain threshold intensity. We used 
a threshold of 4 times of rms sky noise on the region. Our routines count the number of stars 
in each cell of 2 ' x 2' defined on the DSS image and obtain the visual extinction according to 
Dickmanl Jl978h : 



A) 



[l/b)log{no/ni) 



where no is the average number of stars in the comparison cell and n\ is the to- 
tal star number obtained by the count. We obtained star count values in regions away 
from the filamentary structure toward the north-east arou nd (a,5)?nr) n = (0 7 h 44 s ,-41 o 27') 
and computed an average no = 62 for a box of 20'x20 / . Ivan Rhijnl (119291 ) tables yields 
the b parameter and we used b = 0.44. This yalue is the best compromise from aver 



aging several lines of sight tow ard darks clouds (jGregorio-Hetem. Sanzovo fc Lepinelll988 



Andreazza fc Vilas-Boad 119961 ). The maximum error introduced by an incorrect b value is 
not more of 0.3 mag for Ay < 1.5 mag (see below). 

We transformed the extinction obtained from digitized photographic plates (A phc<t ) to 
visual extinction Ay (A5500A) using the normalize extinction curve from iBless fc Savage 
(119721 ). The used DSS image corresponds to Schmidt plates with the emulsion/filter com- 



bination IIIaJ/GC395 centered in 4500A ( Lasker et al.lll990 ). We obtained the necessary 
corre ction by a polynomia l fit in the range 0.4— 3.4/im of the normalized extinction curve 
from lBless fc Savage! (119721 ) and then recalculated the extinction into the photographic band. 
The resulting relation was: 



A, 



0.79A phot 



1 IRAF is distributed by the National Optical Astronomy Observatory which is operated by the Associa- 
tion of Universities for Research in Astronomy Inc., under cooperative agreement with the National Science 
Foundation. 
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We then built an extinction image of 30x 60 cells (== l°x 2°) for the region. The error 
in each cell grows with the extinction value (jDickman! Il978l ). For the cells with no star 
we obtained a lower limit fitting the residuals between consecutive steps of extinction and 
extrapolated the correction for the case n\ = 0. This procedure resulted in assuming n\ 
= 0.28 for cells with no star and the corresponding lower limit of detectable maximum 
extinction was 4.2 mag. Nevertheless, in the regions of interest (i.e. where polarimetric data 
exist) the extinction is relatively low with maximum values of 1.5 mag (^4vmax), well below 
than that limit. Figure [1] shows the extinction maps obtained with this technique. The 
contours go from 0.3 to 1.2 mag in steps of 0.1 mag. 

We must note that procedure used here to estimate the extinction is differential in origin. 
It means that we computed the extiction with respect to the comparison cell uq. In order 
to check if a zero-point bias is p r esent in our extinction maps, we used the reddening maps 
of ISchlegel. Finkbeiner fc David (119981 ). The average extinction in a radius of 5' toward 
the region of the comparison cell was 1.06±0.o43 mag. We consider that this level may 
represent the foreground extinction to the I-front. Therefore, we purposely did not correct 
our extinction maps by this level in order to include just the extinction associated with the 
I-front . 



2 see the Dust Extinction Service available at 



http://irsa.ipac.caltech.edu/applications/DUST/ 



Table 2. Polarization analysis. 



field/7Vcat a 


trend 


Q gauss 


&dgauss 


N 


(P) 


(') 


(8) 


to) 


A6» 


(A9) 2 




(n(H„)) 


B 


A 






(deg) 


(deg) 




(%) 


(%) 


(deg) 


(deg) 


(deg) 


io- 2 


(mag) 


(cm- 3 ) 


(mG) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 



01/14 


1 






14 


1.105 


0.005 


32.9 


02/71 


1 


63.7 


18.1 


66 


0.665 


0.001 


61.2 


03/68 


1 


69.8 


15.2 


43 


0.703 


0.005 


57.9 


04/59 


1 


96.5 


13.2 


35 


1.798 


0.003 


104.0 




2 


40.9 


12.3 


22 


0.870 


0.001 


41.8 


05/76 


1 


40.2 


5.4 


74 


2.284 


0.001 


42.4 


06/68 


1 


10.8 


14.9 


66 


1.723 


0.001 


10.0 


07/61 


1 


171.6 


10.2 


51 


2.529 


0.003 


165.5 




2 


23.5 


2.9 


8 


2.344 


0.012 


23.2 


08/68 


1 


32.5 


8.0 


65 


1.701 


0.003 


32.9 


09/18 


1 






18 


0.429 


0.005 


95.9 


10/78 


1 






78 


0.057 


0.001 


47.3 


11/55 


1 


105.3 


9.5 


43 


0.598 


0.002 


109.4 


12/62 


1 


15.9 


14.3 


29 


1.623 


0.001 


16.0 




2 


116.3 


19.0 


33 


0.799 


0.001 


116.1 


13/50 


1 


104.6 


10.1 


44 


3.113 


0.002 


104.2 


14/63 


1 


75.6 


8.4 


23 


0.435 


0.001 


75.6 




2 


8.1 


19.3 


31 


0.220 


0.001 


9.4 




3 


143.1 


8.9 


9 


0.733 


0.004 


133.0 


15/12 


1 






12 


0.218 


0.003 


158.4 


16/52 


1 






52 


0.251 


0.001 


125.1 



1.2 






0.34±0.24 


280 






1.1 


18.0 


9.9 


0.31±0.10 


255 


0.008 


0.27 


1.8 


15.1 


6.9 


0.23±0.12 


189 


0.009 


0.19 


0.7 


13.2 


5.3 


0.65±0.29 


535 


0.016 


0.28 


1.0 


12.3 


4.6 


0.36±0.17 


296 


0.013 


0.20 


1.0 


5.3 


0.9 


0.42±0.18 


346 


0.033 


0.09 


0.9 


14.9 


6.8 


0.41±0.14 


337 


0.012 


0.25 


0.8 


10.2 


3.1 


0.38±0.10 


313 


0.016 


0.17 


0.7 


2.8 


0.2 


0.45±0.16 


370 


0.064 


0.05 


1.0 


7.9 


1.9 


0.32±0.08 


263 


0.019 


0.12 


1.7 






0.21±0.10 


173 






1.1 






0.18±0.09 


148 






1.2 


9.5 


2.7 


0.33±0.12 


271 


0.016 


0.15 


1.0 


14.3 


6.2 


0.44±0.34 


362 


0.012 


0.25 


0.9 


19.0 


11.0 


0.47±0.31 


387 


0.010 


0.35 


1.0 


10.1 


3.1 


0.42±0.16 


346 


0.017 


0.17 


0.9 


8.3 


2.1 


0.27±0.09 


222 


0.017 


0.12 


1.3 


19.3 


11.3 


0.28±0.14 


230 


0.007 


0.27 


0.8 


8.9 


2.4 


0.36±0.15 


296 


0.018 


0.14 


0.9 






0.58±0.34 


477 






1.1 






0.51±0.12 


420 







a original number of objects in each field from catalog (Paper I). 
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OS 










Fig. 2. — Polarization map for IVS region (see Paper I) with histograms for each position 
analyzed (in numbers). The gaussian fits used to obtain the dispersion of polarization angles 
(AO) in each case are also shown. 
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2.2. Analysis of extinction and masses. 



As noted by IChurchwell et al.1 (119961 ). the IVS far-infrared dust emission follows the 
morphology of the I- front, implying that the radiation field of Vela OB2 and shocks induced 
by the expansion of the shell must heat the dust. These authors also found that the CS 
emission accompanies the I-front and the emission from dust. 

Analyzing Figure [1] we observe that the extinction structures accompany the I-front 
pattern and a good correlation with the extended emission of dust exists (see Figure 2 in 
Paper I). Therefore, it is evident that dust, which originates the observed extinction, was 
shifted to the west of the I-front but accompanying its eastern pattern including the north 
and south curvatures observed on it. Although the extinction values are relatively small (up 
to ~1.5 mag) , the morphology o bserv ed on the extinction map suggests a clumpy region also 
invocated by IChurchwell et al.1 (119961 ) using CS maps. 



In order to better visualize the extinction structures, we choose regions with Ay > 1 
mag. This level is a good compromise (~| ^4vmax) to detect significant structures in our 
maps. Figure [1] shows these structures. We analyse some physical properties of them in 
Table [TJ In columns (2) and (3) we show the equatorial coordinates for the positions with 
the highest extinction observed in each structure. In column (4) is the highest extinction 
value in each structure. Column (5) shows the molecular hydrogen column density estimated 
at the position of highest extinction in e ach structure using the standa rd gas-to-extinction 
ratio (Nn 2 /Ay = 0.94 x 10 21 cm -2 mag~ 1 . lBohlin. Savage fc Drakdl 19781 ). This ratio assumes 
that most of hydrogen is in molecular form and Ry = 3.1. In column (6) is the area in square 
arcminutes associated with each structure and estimated by adding all the cells with Ay > 
1 inside of it. Column (7) shows the dimension associated with the structure in parsecs and 
estimated by -x/ 2 (area). Finally, in column (8) is a lower limit of the mass of each structure 
estimated using iDickmanl (119781 ): 



M = (ad) 2 (N u JAy)mJ2 Ai v, 



where a is the cell size in radians, d is the distance to the cloud in centimeters, (Ah 2 / Ay) 
is the standard gas-to-extinction ratio, m = fin 2 mH = 2.8mjy is the mean particle mass 
(allowing for 10% He by number), /ih 2 is the mean molecular weight with respect to the 
number of H molecules, and A^ is the visual extinction in each reseau element i. 

Among the four structures (A, B, C and D) shown in Figure [H only the A stru cture 
is coincident with the condensation named V12 in IVilas-Boas. Myers, fc Fuller! (119941 . here- 
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i nafte r VMF) in the Vela region. In the southern dark clouds catalogue of lHartley et al. 
( 119861 ). this structure has Ay = 6 mag, and in VMF the extinction is Ay = 3.7 mag using 
C 18 linewidths. Our extinction value (Ay = 1.42 mag) is clearly the lowest. VMF also 
estimated the scale size, hydrogen column density and mass for V12 and obtained 0.15 pc, 
3.4 x 10 21 cm -2 and 1.3M© (lower limit), respectively. The differences with our values in 
Table [1] can be traced to a lower distance used by VMF (300 pc versus 450 pc in this work) 
and the conversion factor of C 18 colu mn density to visua l extinction used by VMF (Ay 
= 6.4 x 10~ 15 iV(C 18 O) + 3.2 mag from lNozawa et al.lll99ll ). In contrast, we obtained our 
extinction values directly from the DSS image. 

The A and B structures in our extinction ma p are coincident w i th A and B clumps seen 
in the distribution study of CS column density of IChurchwell et al.l (119961 ) done in the same 
region. In summary, from Tabled] we estimate a range for typical masses of 1.5—6 M Q and 
a mean length scale L ~0.47 pc for the extinction structures found by us. 
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Fig. 3. — Polarization map overlaying the extinction map. The polarization scale is shown 
on the upper right. The contours are the same as in Fig. [TJ 
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3. Polarization analysis 



In order to better quantify the observed polarization pattern in each one of sixteen fields 
studied in Paper I, we present a statistical analysis in Table [2j Each field is given in column 
(1), together with the catalog stars within the field (iV cat ). We construct polarization angle 
histograms (see Figure [2]) for each field and used gaussian fits to obtain the representative 
polarization angle 9 gauss and its dispersion degauss, indicated in colums (3) and (4). In general, 
one trend prevails in each field but, in some fields, two or three trends appear, and this is 
indicated in column (2). In a few of fields, the star sample was poor, the fit was not possible 
and there is no evident trend. It is interesting to note that none one of the polarization 
trends is parallel to the mean direction of the Galactic Plane {9 ~151°) toward this line of 
sight (/ = 256°, b = —9.2°). This suggests that the polarization data is indeed sampling the 
local magnetic field associated with the IVS region. 

To improve on the precision of the mean polarization values for each trend, we included 
a filter in each field. We select those objects with polarization angle between {6 gauss — 
Zxcegauss, 6 gauss + 2 x ao gauss ) . The number (iV) of stars with P/ap > 10 in each filtered 
subsample is given in column (5). We then estimated the mean Stokes parameters, (Q) and 
(U), for each subsample, from the individuals values for each star Ui) weighted by the 
error (<jj) according to: 



The estimated mean polarization value (P), its associated error (a) and its mean po- 
larization position angle (8) are then given by: 



These are presented in columns (6), (7) and (8). We show the mean error of the polariza- 
tion angle in column (9). This was estimated from (ag) = &8i/N where og i = 28?65 (J Pi /pi 



(Q) = 5>M 2 )/5> r2 



(P) 
(*) 

(0) 



VW+W 2 
(i/E-r 2 ) - 5 



0.5tan-\(U)/(Q)). 
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is the individual error of the polarization position angle (jSerkowskil Il974l ) for an star with 
polarization (pi ± cr Pi ). Column (10) shows the dispersion of polarization angles corrected 



in quadrature by its mean error (A8 
this dispersion (in radians). 



a 



gauss 



{(To) 2 )] column (11) shows the squared of 
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5 % 




RA(2000) 



Fig. 4. — Polarization map overlapping the extinction map for fields 05, 06, 07 and 08. The 
polarization scale is shown up to the right. The contours are the same as in Fig. [TJ 
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40 20 7 43 00 40 20 7 42 00 

RA (2000) 



Fig. 5. — Polarization map overlapping the extinction map for fields 01, 04, 11, 12 and 13. 
The polarization scale is shown up to the right. The contours are the same as in Fig. [TJ The 
condensations A, B and C are also indicated. 
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4. Polarization and extinction 

To help our analysis of correlations between polarization and extinction, we estimated 
the average extinction for objects in each one of the polarization trends (column (12) of 
Table [2]). The position of each star in the polarimetric catalogue is well determined with a 
typical precision of less than 1". We then identified its position in the extinction maps and 
attributed the value of the extinction cell where the star was located. This procedure may 
underestimate the extinction value assigned to a star because the resolution of 2'x2' of the 
extinction maps may not resolve variations of extinction within one extinction cell. 

Figure [3] plots the extinction map of Figure [T] along with the polarization map (Figure 
5 in Paper I) of the region. The polarization vectors, which are from stars with P jop > 
10 in the catalog, are concentrated in regions with significant extinction along the I-front. 
The extinction is low (Ay < 1.5 mag) but is enough to produce the observed polarization. 
The range of extinction sampled by our polarization data is similar to the extinction found 



in the outer layers of darks clouds and in th e diffuse ISM (^4 V < 3 mag, IVrba fc Rydgren 



1985c iKenyon. Dobrzycka fc Hartmannlll994j ) and belo w the threshold extinction of the tran- 



sition from bare to H ^O-mantled grains (Ay ~3.3 mag. iGerakines. Whittet fc Lazarianlll995 



Whittet et al.ll2Q0lh . 



In order to quantify the interval of extinction sampled by our catalogue, we analysed 
two subregions (shown in Figures H] and [5]) with well defined polarization patterns. 

Figure H] shows the four fields (05, 06, 07 and 08) at the top of northern ridge (see 
Figure[2]) around (a, 5)2000 = (07 h 39™5,-41°37 / ). This region presents an observed polarization 
higher than average and a very smooth overall pattern as noted in Paper I. The range of 
mean polarization goes from 1.7% to 2.5% and the mean extinction from 0.32 to 0.45 mag. 
It is interesting to note that the fields 06 and 07 (trend 1), closer to I-front wall and aligned 
with it, present a dispersion of polarization angle (A#) at least twice than ones of the fields 
05 and 08. It reinforces the suggestion of a magnetic field enhancement in the compressed 
postshock g indicated in Paper I. 

The region shown in Figure [5] (fields 01, 04, 11, 12 and 13; see Figure [2]) presents a 
very complex polarization pattern. This complexity seems associated with the A, B and C 
extinction structures. This region has a range of mean polarization of 0.6%— 3.1% with the 
highest value in field 13. The mean extinction goes from 0.33 to 0.47 mag. As we will see 
in £j6j the larger range of polarization values is a indication of significant variations in the 
polarization efficiency in small angular scales in this subregion. 

In some fields, our detection of trends through 9 gaussian fits fails in cases where a 
good polarimetric sample seems to be present. This is evident in fields 09, 10 and 16 (see 
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Figure[2]). Nevertheless, a non-detection of a polarization trend seems correlated with low 
extinction regions as can be verified by the mean extinction of fields 09 and 10 that are 0.21 
mag and 0.19 mag, respectively. However, field 16 has a high mean extinction (0.51 mag) 
but this value may be affected by the contrast gradient present at the border of the DDS 
image used. 



5. Magnetic fields 

As we noted in JD, CF gives an estimate of the magnetic field component (B) on the 
plane of the sky: 



B = (4irp) 1/2 v/A6 



X 



where, B is in gauss (G), p is the density (g cm -3 ), v is the turbulent motion velocity (cm 
s -1 ) and A9 is in radians. In terms of energy densities, the ratio of the kinetic energy density 
(Pkin oc pv 2 ) to magnetic energy density (p ma g oc B 2 ) is the mean square fluctuation of B 
f|Zweibellll990h . so 



{pkin/ P mag 

) oc (A9) . 

Therefore, pkin/Pmag and/or B can be obtained from A8 (column (10) in Table [2]). In 
some fields (01, 09, 10, 15, and 16), it was not possible to obtain information of AO because 
of a poor polarimetric sample or the polarimetric pattern is random. In the remaining 
fields, we estimated (A6) 2 for each trend (column (11) of Table |2j). We can note that (A6) 2 
shows a range of (0.2— 11.3) x 10 -2 . This can be interpreted as variations of (at least) one 
order of magnitude in pkin/p-mag along the region. Nevertheless, this large variation range 
observed in (A6) 2 is consistent with the fact that some regions present random position 
angles. Extreme cases would be fields 09, 10 and 16, where the impossibility of estimating 
A9 may be interpreted as the kinetic energy density prevailing over the magnetic energy 
density. 
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Fig. 6. — (a) Polarization versus vis ual extinction. The solid line represen ts the higher limit 
for ISM polarization [P max = 3 Ay; ISerkowski. Mathewson fc Fordlll975l). The dashed line 
represents the higher limit towards Cha I (P max = 4.5 Ay; IWhittet et al.lll994l ). The dot- 
dashed line is the upper li mit (P = 14 Ay ) for dust grains consisting of completely aligned 
infinite dielectric cylinders ( lWhittetlll992l ) that covers 96.3% of the sample, (b) Polarizing 
efficiency versus extinction. The dashed line represents a least-squares power-law fit to our 



sample (Py/Ay = 2.54±0.08 Ay 



-0.61±0.04 



mag 1 ). 
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Fig. 7. — (a) Polarization versus visual extinction for fields in Figure HI The solid, 
dashed and dot-dashed lines as in Figure [61 (b) Polarizing efficiency versus extinc- 
tion. The dashed line represents a least-squares power-law fit to this sub-sample 
(Py/A Y = 3.50±0.25 ,4 v -o-60±o.o7 % mag -i). 
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Fig. 8. — (a) Polarization versus visual extinction for fields in Figure [5j The solid, 
dashed and dot-dashed lines as in Figure [61 (b) Polarizing efficiency versus extinc- 
tion. The dashed line represents a least-squares power-law fit to this sub-sample 
(Py/A Y = 2.16±0.14 ^ v -o.79±o.or % mag -i). 



-22 - 



In order to calculate the strength of the magnetic field, we assume that the turbulent 
velocity (u) is represented by the dispersion velocity of a Gaussian distribution (AV 2 = 81n2 
v 2 ) where AV is obtained from the mean FWHM of an appropriate molecular line. This 
dispersion may represent the gas, of density p, i n which the polarization arises. The CS emis- 
sion in this region yields AVcs = 1-4 km s _1 (IChurchwell et al.lll996l ). As noted by these 
authors, the density, p = mn(H 2 ), needed for CS excitation requires a particle density n(H 2 ) 
~10 4 cm -3 . As the CS emission samples the higher density cores, its density may represent a 
higher limit for the outer layers of dust clouds which our polarization measurements sample. 
Vilas-Boas. Myers. &: Fuller! (120001 ) analised CO emission in the Vela region and found mean 
dispersions of AVbco = 0.7±0.2 km s _1 and AVcisq = 0.6±0.2 km s _1 . As the CO emission 
probably samples the envelope regions of molecular clouds, its dispersion may better repre- 
sent the turbulent motions. Therefore, we believe that a dispersion of AV = 1.0 km s" 1 is 
a good compromise to consider as input parameter in our computation. 

We used the average extinction (column (12) of Table [2]) to infer the mean particle 
density (n(fi2)) in each polarization trend. We take, 



(n(H 2 )> = (A Y )(N H JA Y ) 1 , 



where lis the typical size (in parsecs) of an extinction cell assuming a distance of 450 pc. 
The computed values are indicated in column (13) of Table [2J The overall average density 
from the trends is 312±97 cm -3 . This value is consistent with our sampling the low density 
regions, as discussed above, and will be used in the following. 

Therefore, with the above parameters and using eq. [TJ the magnetic field expression 
becomes: 



5(G) 



2.88 x 10" 6 ( 



n(H 2 



312 cm" 3 ' 



AV 



1.0 km s 



-)/A6. 



In this way, we estimated B (column (14) of Table El) for each field a nd fo r each trend 
with a well defined dispersion of polarization angle {A8). lOstriker et al.l (120011 ) noted that 
A9 < 25° yields a good estimate of the plane-of-sky magnetic field strength if a multiplying 
factor of ~0.5 is applied to CF formula. As we can note in Table [21 all the A9 values are 
lower than this limit and the above expression also include this correction. As usual, use 
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of the CF formula for estimating B carries an uncertainty from all discussed here. Apart 
from the above mentioned multiplying factor, the uncertainties in n(H 2 ), AV and A6 imply 
a typical uncertainty of at least 50% in B. 

The average strength of the magnetic field along the I-front is [0.018 ± 0. 013] mG 



and th e full range of variation is ~one order of magnitude (0.007—0.064 mG). As IZweibel 



(119901 ) pointed out, if the line of sight samples M different regions with different magnetic 
orientations, B would tend to be overestimated by y/J7. This averaging effect may be 
important, for example, in the region of complex polarization pattern showed in Figure [5] 
where clumpy structures are present. 



Finally, IChurchwell et al.l (119961 ) concluded that the potential energy of this region is just 
~0.05 of the total kinetic energy (thermal plus turbulent) and the expanding shell pressure 
avoids the dissipation of the region. From the estimates of B along with the observed fact 
that, in general, the magnetic pressure dominates over the turbulent pressure in several parts 
of the shell, we conclude that the magnetic component also contributes in an important way 
to the dynamical balance of this region. 



6. Polarizing efficiency 

We combined our polarization data with the extinction maps obtained by star counts 
to analyze the polarizing efficiency of the dust in the region. Figure [6^ plots the visual 
polarization percentage versus the visual extinction for the 856 objects in the catalogue. 

As we noted above, the extinction range in our sample extends up to ~1.5 mag and the 
maximum polarization is ~6%. The solid line in Figure represents the upper limit for 



optim um polarization efficiency in diffuse ISM (P max = 3Av; ISerkowski. Mathewson fc Ford 



19751 ) and just 30.2% of the sample is under this limit. IWhittet et al.l (119941 ) found a higher 
upper limit for optimum polarization (P max = 4.5Ay) toward some lines of sight in the Cha I 
dark cloud. They concluded that a high degree of alignment is present with the magnetic field 
lines essentially perpendicular to our line of sight. If we consider this upper limit (dashed 
line in Figure [6^) the percentage of our sample under this limit rises to 49.8%. 

It seems to be clear that part of the dust that produces the observed alignment in 
this region has different properties than diffuse ISM dust and it is especially true for low 
extinction (Ay < 0.8 mag) regions. On the other hand, the higher polarization efficiency 
for low extinction can be a indicator that an optimum alignment is present in some regions 
of I-front maybe favored by a privileged geometric view (i.e. parallel to plane of the sky) 
of the magnetic field with respect to the wall seen edge-on. This fact can be quantified by 
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the theoretical upper limit (Py/Ay < 14 ) for dust grains consisting of completely aligned 



infinite dielectric cylinders ( )Whittetlll992l ) that covers 96.3% of the sample (dot-dashed line 
in Figure Eh). 

Nevertheless, and as mentioned in §HJ the possibility of our 2'x2' extinction cell will be 
underestimating the extinction can not be discarded. If this is the case, higher extinction 
values can be masked by the observed clumpy structures and an optimum alignment better 
than different dust properties is more plausible to explain the correlations found. 
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logW(H 2 ) 



Fig. 9. — Observed mass-to-magnetic flux ratio in units of the critical value (corrected by 
bias projection) versus molecular column density. In black dots are the As evaluated in each 
polarization trend (with B estimated) from Table [2j The IVS region is dominated by the 
subcriti cal regime (A < 1). The dashed region represents all the As available for molecular 
clouds (jChrutcherl 120041 1 . The box inside the figure is a zoom to show more detail. 
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In Figure (6jo, we plot the polarizing efficiency versus the extinction and it is very clear 
as regions with low extinction are more efficient to polarize the light of background stars 
than hi gh extinction regions. This fact already has been obse r ved previously in darks 



1992 



1995 



clouds (IVrba. Strom, Stroml 1 19761: IVrba. Coyne. &: Tapia 



Vrba. Covne. fc Tapialll993l : McGregor et al.lll994j : IWhittet et al 



1981 



Vr 



ja. Marraco k, Strom 



1994 



HH^^HH^H^HHHdv" 1 ■ O ' — " « ^ 7 '— — " « ^ 7 ■ " — 

Goodman et al.lll995l ). However, in I- fronts of shells in expansion, as this study, pre 



Gerakines. Whittet fc Lazarh 



vious evidence does not exist in literature. 

The dashed line in Figure [6b shows a least-squares power-law fit to our sample (Py/Ay 
= 2.54±0.08 Ay~°' 61±0M % mag -1 ). The depende nce of the p olarizing efficiency with the 
extinction found by us is very close to the models of lJoned (119891 . Pj A oc A~ - 5 ) that assumes 
a magnetic field with random and unif orm components. Our fit also compares very well with 
Gerakines. Whittet fc Lazarian Jl995l . Pj A oc A a56±0 - 17 ) for observed data in Taurus dark 
cloud. 



The ratio Pj A may be written as Pj A oc A 2k 1 (IGerakines. Whittet fc Lazarianlll995 ). 
wher e k represents the amplification of the magnetic field with a gas density (jMouschovias 



1978 



B oc n k ). Our fit y i elds k = 0.195. This value is below the observed lower limit 
(k = 0. 33) of iMouschoviasi (119781) and the mode ls of magnetic fields in turbulen t mole cular 
clouds (jOstriker et al.ll200ll ). As pointed out by IGerakines. Whittet &: Lazarian! (119951 ) this 
discrepancy could be result of ignoring effects of small-scales inhomogenties in the magnetic 
field and of the coupling of gas and dust tem peratures in regions of high dens ity in the 



assumed polarizing efficiency [Pj A oc P> 2 / n from lVrba. Covne. fc Tar>ialll98ll. Il99sh. F inally , 
the coefficient of our fit (2.54) is close to Pj A = 3 (from lSerkowski. Mathewson fc Fordlll975l ) 
and this suggests that the higher extinction regions (0.8 < Ay < 1.5) of our sample are more 
normal and comparable to diffuse ISM. 

As comparison, we made the same as analysis for the two subsamples shown in Figures H] 
and The results are shown in Figures [7] and [HI respectively. As we can note, the region 
sampled in Figure H] corresponds to low extinctions, below ~0.8 mag. We find a high polar- 
izing efficiency and this subsample seems to be representative of the overall trend observed 
in the full sample (Figure ED • Note that the exponents of the fits for this subsample and 
the full sample are essentially the same. On the other hand, the very complex polarization 
pattern shown in the region indicated in Figure [5] has a higher extinction range (up to ~1.5 
mag, Figure [8]) and the polarizing efficiency is lower if compared with the previous region. 
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Mass-to-flux ratio 



Cardelli Savage! (119881 ) noted that IRAS point sources are coincident with the densest 
portion of this region. They speculated that shocks might have triggered local collapse and 
possibly star formation. However, they associated the clumpy structure seen in this region 
with clumps of warm dust with peak emission in 60/zm band and star formation was ruled 
out. 

The importance of magnetic fields in the evolution of interstellar clouds and star for- 
mation can be tested if the mass-to-magnet ic flux ratio can be qu antified (M/$). We can 
evaluate this issue using the A parameter (jChrutcher et al.l I2004J . and references therein) 
defined as: 



A = 



(M/$) 
(M/$) 



actual 



7.6 x 10 



crit 



-2lgg2) 

B 



(2) 



where (M/Q) crit is the critical value for the mass that can be supported by a magnetic 
flux, 7V(H 2 ) is in cm -2 and B is in fiG. In the subcritical regime (A < 1), magnetic support 
prevents the collapse and the star formation and probably this is the case in the IVS region. 

We evaluated the A parameter (using eq. [2]) for each polarization trend where B has been 
estimated. The results are indicated in the last column of Table |2J We obtained iV(H 2 ) from 
the mean extinction (Ay) for each trend (column (12) in Table [2]) and using the standard 
gas-to-extinction ratio (as in §2.21) . As (Ay) and B are estimated in exactly the same area, 
this guarantees a proper determination of A. Clearly, all the regions with an alignment trend 
are in the subcritical regime (A < 1) between 0.05 and 0.35. If we ap plied the 1/3 correction 
factor to take in count the projection bias (jChrutcher et al.l 12004 ). the subcritical regime 
is even sharper. This is evident in Figure [9] that shows the distribution of A (corrected by 
projection bias) as a function of the molecular column density. 

It seems clear that magnetic support is important in the IVS region and the data are 
consistent with the prediction of ambipolar diffusion models with cloud envelopes being ini- 
tia lly subcritical. Thi s is als o consistent with non-star formation in this region as speculated 
bv lCardelli fc Savage! Jl988h . 



It is interesting to note that the two points with the largest column densities and largest 
As (Figure [H]) corresp ond to polar i zation trends associated with the high extinction structures 
A and B (see §2.2p . IChrutcherl (120041 ) analyzed all the As available for molecular clouds 
(dashed region in Figure [9]) and noted a slight indication that for large colum densities, A 
may be supercritical, and for small column densities, subcritical. Alhough that our data in 
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Figure [9] samples just a interval of lower iV(H2) values (~10 20 cm 2 ), evidence of A rises with 
iV(H 2 ) can not be discarded. 

8. Conclusions 

Our conclusions can be summarized as follows: 

1. Extinction maps obtained using automatic star counts toward the western side of 
IRAS Vela Shell yield evidence of clumpy structures with typical lower limits of masses 
between 1.5 and 6 M and a length scale L ~0.47 pc (assumed a distance of 450 pc). 

2. We estimated the strength of the local magnetic field on the plane of the sky using the 
dispersion of polarization angles and educated estimates for other cloud properties. The mean 
value computed along the I-front was [0.018±0.013] mG and the full range of variation was 
(0.007—0.064) mG. In terms of energy densities, the magnetic pressure generally dominates 
over the turbulent motions along the I-front. In a few cases, the polarization angle appears 
to change randomly, suggesting that the kinetic energy can dominate. 

3. We also investigated the polarizing efficiency combining the polarization data with 
the extinction maps. We found high polarizing efficiency in low extinction regions. It can be 
explained by different properties of dust as compared with the general ISM dust. Another 
more plausible possibility to explain this high polarizing efficiency is an optimum mechanism 
of dust alignment maybe favored by the magnetic field lines are nearly parallel to the plane 
of the sky and perpendicular to the I-front. 

4. Using the strength of magnetic field and the extinction estimated in each polarization 
trend was possible to evaluate the mass-to-magnetic flux ratio along the IVS region. The 
range of the A parameter we found (0.05—0.35) confirms that the magnetic support is domi- 
nant along this I-front and a subcritical regime prevails. The A— N(H 2 ) relation we find joins 
smoothly the overall relation previously obtained for higher density regions. This provides 
general support for the evolution of initially subcritical clouds to an eventual supercritical 
stage. 

The authors wish to thank the anonymous referee for his/her careful reading and sug- 
gestions that helped to improve the paper. AP is thankful to CAPES and FAPESP (grant 
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